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Abstract: Detecting near-surface soil freeze-thaw cycles in high-altitude cold regions is important for 
understanding the Earth’s surface system, but such studies are rare. In this study, we detected the 
spatial-temporal variations in near-surface soil freeze-thaw cycles in the source region of the Yellow River 
(SRYR) during the period 2002—2011 based on data from the Advanced Microwave Scanning Radiometer 
for the Earth Observing System (AMSR-E). Moreover, the trends of onset dates and durations of the soil 
freeze-thaw cycles under different stages were also analyzed. Results showed that the thresholds of 
daytime and nighttime brightness temperatures of the freeze-thaw algorithm for the SRYR were 257.59 
and 261.28 K, respectively. At the spatial scale, the daily frozen surface (DFS) area and the daily surface 
freeze-thaw cycle surface (DFTS) area decreased by 0.08% and 0.25%, respectively, and the daily thawed 
surface (DTS) area increased by 0.36%. At the temporal scale, the dates of the onset of thawing and 
complete thawing advanced by 3.10 (+1.4) and 2.46 (+1.4) days, respectively; and the dates of the onset 
of freezing and complete freezing were delayed by 0.9 (+1.4) and 1.6 (£1.1) days, respectively. The 
duration of thawing increased by 0.72 (+0.21) day/a and the duration of freezing decreased by 0.52 
(£0.26) day/a. In conclusion, increases in the annual minimum temperature and winter air temperature are 
the main factors for the advanced thawing and delayed freezing and for the increase in the duration of 
thawing and the decrease in the duration of freezing in the SRYR. 
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1 Introduction 


The near-surface soil freeze-thaw cycle is related to the fluctuations in soil temperature around 
0°C due to seasonal or diurnal temperature changes. Such fluctuations are common at high 
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latitudes or within high-altitude zones. Approximately 50x10° km?, or 55% of the land surface 
soil in the Northern Hemisphere, undergoes an annual soil freeze-thaw cycle (Zhang et al., 2009). 
The mechanism of the soil freeze-thaw cycle involves alternating water-ice phase transitions in 
soil; this phase transitions act as “switch” in nature. That is, when soil is thawed, the “switch” 
opens, and hydrological, ecological and other land surface processes become active; when soil 
freezes, the “switch” closes, and the land surface processes enter a period of dormancy or slow 
activity (Zhang et al., 2011; Yi et al., 2014). Changes in the near-surface soil freeze-thaw cycle 
strongly affect the climate system, the hydrological cycle, the carbon and nitrogen cycles, surface 
energy fluxes, plant growth, soil erosion, and other land surface processes (Poutou et al., 2004; 
Wang et al., 2006; Henry et al., 2008; Matzner and Borken, 2008; Sun et al., 2008; Nishimura et 
al., 2011; Buermann et al., 2013; Wang et al., 2015). 

Previous research on the monitoring of near-surface soil freeze-thaw cycles has mainly been 
carried out using data measured at individual sites (Yang et al., 2007; Li et al., 2012; Guo and 
Wang, 2014). However, this traditional method is time consuming and cannot meet the 
requirements of large-scale and dynamic monitoring because the measurement sites are often 
sparse and unrepresentative. Satellite-based remote sensing provides an efficient mean to detect 
near-surface soil freeze-thaw status over long periods of time at global, continental, and regional 
scales (Zhang et al., 2003; Bartsch et al., 2007). Especially, passive microwave remote sensing, 
which has a high temporal resolution and all-weather capabilities, is suitable for this task (Zhang 
et al., 2011; Zhao et al., 2011; Li et al., 2012; Singh et al., 2013). The main sources of passive 
microwave remote sensing data include the Scanning Multichannel Microwave Radiometer 
(SMMR, 1978-1987), the Special Sensor Microwave Imager (SMM/I, 1987-present) and the 
Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E, 
2002-2011). Compared with SMMR and SMM/I, AMSR-E is newer and more efficient ifor 
detecting soil freeze-thaw status, although relatively few studies have used the data collected by 
this sensor. 

Researchers have developed several methods to monitor near-surface soil freeze-thaw status in 
the application of passive microwave remote sensing data, i.e., the dual-indices algorithm, change 
detection algorithm (Smith et al., 2004), decision tree algorithm and discriminant function 
algorithm (Zhang et al., 2003; Li et al., 2012). It is worthy noting that the dual-indices algorithm 
has been widely applied to study the freeze-thaw patterns (Han et al., 2010; Xie et al., 2013). This 
algorithm uses a 37-GHz brightness temperature and a negative spectral gradient between the 18 
and 37 GHz brightness temperatures to detect soil freeze-thaw status (Zuerndorfer and England, 
1992; Judge et al., 1997). Moreover, the thresholds of the dual-indices algorithm must be 
separately determined for each region using data collected from different sites and verified 
through the statistical analysis (Chai et al., 2014). 

The source region of the Yellow River (SRYR), which lies in the interior of the Qinghai-Tibet 
Plateau (QTP) and is known as a “sensitive area” and an “amplifier” for global warming, has 
become warmer at a fast rate than other places on the plateau (Lan et al., 2016). Many lines of 
evidence indicate that climate warming has resulted in a serious degradation of frozen ground in 
recent years over the SRYR, leading to reduction in permafrost area (Wang et al., 2001; Jin et al., 
2009; Matsuoka et al., 2009; Luo et al., 2016). However, the near-surface soil freeze-thaw cycle 
over the SRYR has not been systematically investigated. Therefore, we investigated the 
near-surface soil freeze-thaw cycle in the SRYR during the period 2002-2011 based on the 
AMSR-E passive microwave remote sensing data. The aims of this study were (1) to obtain the 
thresholds of the dual-indices algorithm for soil freeze-thaw state detection in the SRYR; (2) to 
investigate the spatial distribution patterns of the near-surface soil freeze-thaw state; and (3) to 
analyze the variations in the onset dates and durations of the near-surface soil freeze-thaw cycle 
during the period 2002-2011. 


JOURNAL OF ARID LAND 


2 Materials and methods 


2.1 Study area 


The SRYR catchment area, situated in the upper stream of the Longyangxia Reservoir, covers 
approximately 145.3x10° km? in area. Geographically, it is located between 32°10'—36°59’N and 
95°54'—103°24'E, and the average elevation is greater than 4000 m (Fig. 1). The study area is 
characterized by a typical continental plateau climate, with an annual average air temperature of 
—3.0°C to —4.1°C and an average annual precipitation of 300-700 mm. Three landform types exist 
in the SRYR: mountains, hills, and plains. Large areas in the eastern and central SRYR are 
covered with seasonally frozen soil interspersed with isolated patches of sporadically continuous 
permafrost. The permafrost occurs primarily on high mountains in the western part of the region 
(Wen et al., 2011). 
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Fig. 1 Topographic map of the source region of the Yellow River (SRYR) 


We divided the study area into five regions according to the Qinghai-Tibetan Plateau natural 
zonation scheme (Zheng et al., 1979) for analyzing the spatial characteristics of the near-surface 
soil freeze-thaw cycle. These regions are the wide valley basin of the Yellow River source region 
(I), the Zoige plateau (II), the Golog and Yushu plateau gullies (III), the Huangnan Mountains 
(IV), and the eastern margin of the Qaidam Mountains (V) (Fig. 2). 


2.2 Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) 
data 


We obtained the daily AMSR-E passive microwave remote sensing data from the National Snow 
and Ice Data Center (NSIDC) at the University of Colorado, USA. These data are provided in the 
Equal-Area Scalable Earth Grid projection (EASE-Grid) and have a spatial resolution of 25 km 
for all channels. The data are available twice daily from an ascending track (the satellite passed 
over the study region at 13:30 local time) and a descending track (the satellite passed over the 
study region at 01:30 local time) for 10 consecutive years from 19 June, 2002 to 3 October, 2011. 
The radiometer channel frequencies are 6.9, 10.7, 18.7, 23.8, 36.5 and 89.0 GHz, respectively, 
and each frequency is associated with dual vertically and horizontally polarized radiation 
measurements. In this study, the dual-indices algorithm was applied to the 18.5 and 36.5 GHz 
vertically polarized brightness temperature data (Tz_18.sv and Tg_36.sv) from both the ascending and 
descending tracks. 

Some of the AMSR-E data were missing because of the unstable states of the radiometer and 
atmospheric conditions. The missing data included the 211’-217" days and the 256'"-262"™ days 
in 2002, the 304'"-308" days in 2003, the 324" day in 2004 and the 332™ day in 2007. The daily 
AMSR-E data did not cover the entire study area at the middle and low latitudes; therefore, we 
had to splice the data from two adjacent days into one to cover the entire SRYR. 
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Fig. 2 Land cover and natural zonation of the SRYR. The purple solid lines represent the natural zonation 
boundaries. Natural zones: wide valley basin of the Yellow River source region (I), Zoige plateau (II), Golog and 
Yushu plateau gullies (III), Huangnan Mountains (IV), and eastern margin of the Qaidam Mountains (V). The red 
dotted line A represents the north-south profile line at 99.8°E for Figure 7, and the red dotted line B represents the 
east-west profile line at 34.5°N for Figure 6. 


2.3 Meteorological data 


In this study, a land surface temperature (LST) dataset that consisted of data from eleven 
meteorological stations in the SRYR obtained from the Chinese Meteorological Science 
Data-sharing Service (Fig. 1) was used to calibrate and validate the thresholds of the dual-indices 
algorithm. The LST dataset included the daily average land surface temperature (LSTave), the 
maximum land surface temperature (LSTmax) and the minimum land surface temperature (LST min) 
for the period 2002-2011. In addition to the LST data, the monthly average air temperature, the 
maximum air temperature, and the minimum air temperature were also collected for the period 
2002-2011. Meteorological data from five meteorological stations, i.e., Qabqa, Guinan, Golog, 
Zoige and Madoi stations, were used to calibrate the dual-indices algorithm, and data from other 
stations (Xinghai, Henan, Maqu, Hongyuan, Jigzhi and Darlag stations) were used to validate the 
algorithm. 


2.4 Digital Elevation Model (DEM) and land cover data 


In this study, the Advanced Spaceborne Theemal Emission and Reflection Radiometer Global 
Digital Elevation Model (ASTER GDEM), which has a spatial resolution of 30 m, was obtained 
from the Geospatial Data Cloud. A land cover map for 2009, which has a spatial resolution of 300 
m and uses the UN Food and Agriculture Organization’s land cover classification scheme, was 
obtained from the Cold and Arid Regions Sciences Data Center at Lanzhou 
(http://westdc.westgis.ac.cn). The above datasets were pre-processed using a system radiation 
correction and ground control point correction. 

Profile lines can be used to analyze the spatial distribution of onset dates and durations of 
freeze-thaw cycles at regional scales (Bartsch et al., 2007). We created an east-west profile line at 
34.5°N that crossed areas I, II, and III and a north-south profile line at 99.8°E that crossed areas 
IV and III. These lines intersected at the center of the SRYR and passed through areas with typical 
and varied landscapes based on the land use and natural zonation in the SRYR (Fig. 2). 
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2.5 Analysis method 


We used the trend line method to analyze the changes in onset date and duration of the 
freeze-thaw cycles estimated using map algebra in ArcGIS. The computational formula is as 
follows: 


nx xx) OIL) 


slope — n n 
nx -È 
j=l j=l 


where n is the number of years in the monitoring period (n=10 in this study); xj, the quantity 
being studied in year j (onset dates or duration of the freeze-thaw cycle); and 0, the slope of the 
trend line (when @0, the quantity being studied increases over the n years, and vice versa). 

After the trend maps were created, the mean values were computed for each pixel, and the 
standard deviations were calculated as the average deviation of each pixel from the mean value. 
These calculations were implemented using the Spatial Analyst Toolbox in ArcGIS. 


(1) 
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3 Results 


3.1 Calibration and validation of the dual-indices algorithm 


The dual-indices algorithm uses the 36.5 GHz vertical-polarization brightness temperature 
(Tg _36.5v), the spectral gradient (SG) between the 36.5 GHz (Tg 36.sv) and 18.5 GHz (Teg_18.5v) 
brightness temperatures and data from AMSR-E. This algorithm can be expressed as: 


Tg 365v < T atott ? (2) 
SG= Ty 365v -Tp igsv <0, (3) 
36.5-18.5 


where Tutor is the threshold determined from the linear regressions between LST and Tg _36.sv for 
determining the near-surface soil freeze-thaw status. 

LSTs tend to be close to the maximum values at the ascending local time (13:30) for AMSR-E 
and approach to the minimum values at the descending local time (01:30) (Xie et al., 2013; Chai 
et al., 2014). Therefore, we used LSTmax and LSTmin to calibrate and verify the classification 
algorithm for the ascending and descending tracks, respectively. Two linear regressions between 
the LSTs and Tz_36.5sv were established for the period from 1 January, 2007 to 31 December, 2009 
(Fig. 3), and robust correlations (R?=0.772 and R?=0.767, respectively) were identified between 
these quantities and the standard deviations were 3.87 and 3.32 K, respectively. According to the 
regression lines, the thresholds of the dual-indices algorithm were 257.59 K during the day and 
261.28 K at night, which corresponded to the real soil freezing point (273.15 K). A pixel of the 
AMSR-E data was assigned to the thawed status during the day if the ascending Tz 36.sy and SG 
exceeded 257.59 and 0 K, respectively; otherwise, the pixel was assumed to be frozen. A pixel 
was classified into the frozen status at night if the descending Tg 36sv and SG were less than 
261.28 and 0 K, respectively; otherwise, the pixel was assumed to be thawed. Additionally, large 
bodies of water (lakes and rivers), areas of permanent snow and ice and deserts influence 
brightness temperatures and could be mistakenly classified. Therefore, areas up to 50 km with 
those land cover types (e.g., lake, river, snow, ice, and desert) were masked out in this study. 

The following variables were used to evaluate the accuracy of the algorithm: (a) Fy represents 
the number of correctly classified daily total frozen pixels; (b) Fa, the actual number of daily total 
frozen pixels; (c) Ty, the number of correctly classified daily total thawed pixels; (d) Ta, the actual 
number of daily total thawed pixels; (e) FTC, the number of correctly classified daily 
freeze-thaw cycle pixels; and (f) FTC,, the actual total number of daily freeze-thaw cycle pixels. 

The calculation method for the evaluation variables was as follows: (1) Fv/Fa ratio can be used 
for detecting the classification accuracy of the daily total frozen surface; (2) T\/T, ratio, the 
classification accuracy of the daily total thawed surface; (3) FTC\/FTC, ratio, the classification 
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accuracy of the daily freeze-thaw cycle surface; and (4) (Fy+T\+FTCy)/(FatTit FTC.) ratio, the 
overall classification accuracy. 
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Fig. 3 Relationship between land surface temperature (LST) and brightness temperature at 36.5 GHz (Ts-36.5v) 
with vertical polarization based on Advanced Microwave Scanning Radiometer for the Earth Observing System 
(AMSR-E) data for daytime (a) and nighttime (b) passes. LSTmax, the maximum land surface temperature; LSTmin, 
the minimum land surface temperature; RMSE, root mean square error; F, degree of freedom. 


We verified the algorithm with LSTmax and LSTmin measurements from six meteorological 
stations (Xinghai, Henan, Maqu, Hongyuan, Jigzhi and Darlag) during the period 2007—2009. 
Additionally, the same thresholds were used in the calibration algorithm. Based on 13,152 
validated data points, the classification accuracies of the daily total frozen surface, the daily total 
thawed surface, and the daily freeze-thaw cycle surface were 80.46%, 83.65%, and 78.07%, 
respectively, and the overall classification accuracy was 80.55%. These results indicate that the 
dual-indices algorithm can be used to detect the near-surface soil freeze-thaw cycle in the SRYR 
with sufficient accuracy. 

Based on the aforementioned dual-indices algorithm and thresholds, we detected the daily 
near-surface freeze-thaw status derived from the AMSR-E data during the day and at night and 
then classified the status into three types. These types corresponded to the daily frozen surface 
(DFS) where the land surface was classified into the frozen status during both day and night; the 
daily thawed surface (DTS) where the land surface was classified into the thawed status during 
both day and night; and the daily surface freeze-thaw cycle (DFTS) where the land surface was 
classified into thawed during the day and frozen at night. Based on the time series information of 
the daily near-surface freeze-thaw status, four soil freeze-thaw stages in a year were defined: a 
spring freeze-thaw transition stage, a complete thawing stage, an autumn freeze-thaw transition 
stage and a complete freezing stage. If the surface belonged to one of three daily freeze-thaw 
statuses for three successive days, the periods of the surface status were determined as one of the 
four stages, and the first day of these three days was recorded as the onset date of the next 
freeze-thaw stage. 


3.2 Spatial distribution of near-surface soil freeze-thaw status 


Figure 4 shows the monthly variations in near-surface soil freeze-thaw status from September 
2010 to August 2011. As the temperature decreased after September, a DFS occurred near 
Ngroing Lake and Gyaring Lake in area I, while a DFTS occurred near Maduo County and Requ 
River. The DFS and DFTS subsequently expanded, moved southeastward, and gradually occupied 
area III. In January, which is the coldest month of the year, the DFS and DFTS were bounded by 
the A’nyémagén Mountains, whereas the DFTS and DTS were bounded by the Xiqing Mountains. 
At that time, area I, which lies to the west of the SRYR, displayed a DFS status; areas III, IV and 
V, which lie in the middle of the SRYR, displayed a DFTS status; and area II, which lies to the 
southeast of the SRYR, displayed a DTS status. As the air temperature gradually increased, the 
western boundary of the DFS, the eastern boundary of the DTS, and the eastern and western 
boundaries of the DFTS moved westward. Later in March, the DFS disappeared, and a DFTS 
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replaced the DFS in area I to the west of the A’nyémaqén Mountains. The DTS gradually 
increased and then occupied the entire SRYR in July and August. Therefore, the SRYR, except for 
the area with DTS throughout the year, underwent freeze-thaw cycles over long periods of time. 
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Fig. 4 Monthly variations in near-surface soil freeze-thaw status from September 2010 to August 2011 in the 
SRYR 


The areas occupied by the three freeze-thaw statuses showed different trends from 2003 to 
2010 (Fig. 5). The DFS gradually decreased from January to June and began to increase in 
September. Additionally, the peak of the curve for the month with the largest DFS moved 
backward from 2003 to 2010 (Fig. 5a), and the average DFS decreased by approximately 0.08%. 
In spring and autumn, the surface temperature dropped below zero at night, whereas the surface 
temperature exceeded zero during the day, resulting in freezing at night and thawing during the 
day. Being similar with the DFS, the DFTS occurred in autumn or winter, although its peak 
occurred during the early spring (Fig. 5b); and the surface gradually decreased and completely 
disappeared in summer. The DTS increased from January to July or August and then gradually 
decreased (Fig. 5c). 
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Fig. 5 Monthly variations in the daily frozen surface (DFS, a), daily freeze-thaw cycle surface (DFTS, b) and 
daily thawed surface (DTS, c) during the period 2003-2010 
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3.3 Changes in the onset date and duration of the freeze-thaw cycles during the period 
2002-2011 


We analyzed the spatial distribution and inter-annual variation of the onset date of the four stages 
and the duration of the freeze-thaw cycles using the east-west profile line at 34.5°N (Fig. 6) and 
the north-south profile line at 99.8°E (Fig. 7). For the east-west profile line at 34.5°N, the dates of 
the onset of thawing (10—-65" days) (Fig. 6a) and complete thawing (28-115 days) (Fig. 6b) 
moved forward from east to west as the elevation decreased, and the dates of the onset of freezing 
(260-320" days) (Fig. 6c) and complete freezing (310-360" days) (Fig. 6d) moved backward. 
The duration of the spring and autumn freeze-thaw transitions, which were determined by the 
onset of thawing/freezing and the onset of complete thawing/freezing, decreased from east to west. 
The duration was shorter in area II (spring: 45 days; autumn: 28 days) than in area III (spring: 52 
days; autumn: 31 days); moreover, this duration was shorter than in area I (spring: 57 days; 
autumn: 35 days). The dates of the onset of thawing and complete thawing tended to occur earlier 
over the study period (2002-2011), and the dates of the onset of freezing and complete freezing 
tended to occur later over the study period. It should be noted that part of area II did not 
completely freeze in winter after 2006. 
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Fig. 6 Changes in the onset date and duration of the freeze-thaw cycles for the east-west profile line at 34.5°N 


during the period 2002-2011. Onset date of thawing (a), date of complete thawing (b), onset date of freezing (c), 
date of complete freezing (d), and elevations in areas I, II, and III (e). 
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For the north-south profile line at 99.8°E, the dates of the onset of thawing (10'—55" days) (Fig. 
7a) and complete thawing (71‘-104" days) (Fig. 7b) moved forward first and then backward from 
north to south. Similarly, the dates of the onset of freezing (280-318 days) (Fig. 7c) and 
complete freezing (310'-345" days) (Fig. 7d) moved backward first and then forward. Turning 
point is at the transition zone between areas IV and III in the rolling mountain region of 
A’nyémagén. The durations of the spring and autumn freeze-thaw transitions, which were 
disrupted by the elevated mountains, first increased and then decreased along this line from north 
to south. In area IV in the northern SRYR, the spring duration (approximately 52 days) and 
autumn duration (approximately 38 days) of freeze-thaw transitions increased from north to south, 
whereas an opposite trend was observed in area III. Similar to the east-west profile, the date of the 
onset of the freeze-thaw cycles along the north-south line showed a trend toward earlier thawing 
and later freezing from 2002 to 2011 against the background of climate warming in the SRYR. 
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during the period 2002-2011. Onset date of thawing (a), date of complete thawing (b), onset date of freezing (c), 
date of complete freezing (d), and elevations in areas III and IV (e). 


Figure 8 shows the spatial distribution of the onset date and duration of the freeze-thaw cycles. 
The dates of onset of thawing advanced or moved forward by 3.1 (+1.4) days over the SRYR (Fig. 
8a), and most of these advanced days appeared in area IV, followed by area III.The complete 
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thawing date (Fig. 8c) was advanced by 2.46 (+1.4) days and most of these advanced days 
appeared in area II, with a decreasing tendency in the number of advanced days from southeast to 
northwest. The dates of onset of freezing (Fig. 8e) and complete freezing (Fig. 8 g) were delayed 
by 0.9 (41.4) and 1.6 (41.1) days, respectively. Most of these delayed days appeared in parts of 
areas II and IV in the eastern SRYR, whereas no significant changes occurred in area I. 
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Fig.8 Spatial distribution of trends in the onset date and duration of the freeze-thaw cycles (in days). Onset date 
of thawing (a), duration of the near-surface soil freeze-thaw transition during spring (b), date of complete thawing 
(c), duration of the near-surface soil thawing (d), Onset date of freezing (e), duration of the near-surface soil 
freeze-thaw transition during autumn (f), date of complete freezing (g), and duration of the near-surface soil 
freezing (h) in the SRYR during the period 2002-2011. 


The duration of the spring freeze-thaw transition (Fig. 8b) in the middle of the SRYR increased 
by 2.1 (+0.91) days, whereas that in western Ngoring Lake, Zoige Wetland and the northernmost 
SRYR the duration decreased by 1.21 (+0.56) days. Similarly, the duration of the autumn 
freeze-thaw transition (Fig. 8f) in the central region increased by 1.1 (+0.42) days, and the 
transition decreased by 1.21 (+0.56) days in the other parts. Overall, the duration of complete 
thawing (Fig. 8d) increased by 0.72 (+0.21) days/a with the maximum increase occurred in area 
IV and the minimum increase in western Ngoring Lake and northernmost area. Overall, the 
duration of complete freezing (Fig. 8h) decreased at a rate of 0.52 (+0.26) days/a with the greatest 
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decrease occurred in the northeastern SRYR at lower altitudes and the smallest increase in the 
southwest SRYR at higher altitude. The above results suggest that the changes in the freeze-thaw 
cycles were most severe in areas III and IV in the northeastern and southeastern parts of the 
SRYR in relatively low-lying regions. 


4 Discussion 


Many researchers have found that an increase in air temperature, especially in the cold season, 
results in the shortening of the freeze duration and the lengthening of the thawing duration (Li et 
al., 2012b; Guo and Wang, 2014; Wang et al., 2015). In this study, the annual average air 
temperature (ATaye) of the SRYR during the period 2002-2011 increased at a rate of 0.023°C/a 
(Fig. 9). Moreover, the increase rate of the annual minimum air temperature (AT min, 0.066°C/a) 
was higher than that of the annual maximum air temperature (ATimax, 0.042°C/a), and the increase 
rate of average air temperature in winter (AT winter, 0.072°C/a) was higher than that of other 
seasons, suggesting that there was an asymmetric increase in air temperature in different times 
and seasons (Zhao et al., 2004; Liu et al., 2006). That is to say, the increase in air temperature was 
higher in cold seasons and at night than in warm seasons and during the day, respectively. As shown 
in Table 1, ATmin was significantly correlated with the durations of complete thawing and complete 
freezing, and ATwiner Was significantly negatively correlated with the durations of the spring 
freeze-thaw transition and complete freezing. These results confirm that the increases in ATmin and 
ATwinter resulted in the advanced thawing and delayed freezing, and the increase in the thawing 
duration and the decrease in the freezing duration. 

With respect to variations in the near-surface soil freeze-thaw cycles from other parts of the 
world, similar results have been found in Mongolia (Han et al., 2010), central USA (Anandhi et 
al., 2013) and China (Wang et al., 2015). However, results from this study indicate that the onsets 
of thawing during the period 2002-2011 in the SRYR advanced by approximately 0.34 day/a, 
whereas it advanced by 0.15 day/a in all of China (Wang et al., 2015) and 0.10 day/a in Kansas, 
USA (Anandhi et al., 2013). Similarly, the duration of freezing exhibited a larger variation in the 
SRYR (approximately 0.52 day/a) than in all of China (approximately 0.25 day/a) (Wang et al., 
2015) and in Kansas, USA (0.13 day/a) (Anandhi et al., 2013), but less than that in Mongolia 
(1.40 day/a) (Han et al., 2010). The differences among these results should not be ignored and 
may be attributable to two major reasons. First, as previously stated, increases in air temperatures, 
particularly the ATmin and ATwinter, significantly affect variations in near-surface soil freeze-thaw 
cycles. It is clear that the changing rates of air temperature differ among study areas. For example, 
the increase rate of air temperature was faster in the SRYR than in all of China (Lan et al., 2016) 
and also in central USA (Brunsell et al., 2010) and slightly slower than in Mongolia (Han et al., 
2010). Second, the use of different detection methods for the near-surface soil freeze-thaw cycle 
may lead to different results. Han et al. (2010) and Li et al. (2012) used soil temperature measured 
at 06:00 to represent the onset or offset of the soil freezing events and temperature measured at 
18:00 to represent the onset or offset of the soil thawing events based on the SMM/I data. 
However, this may have ignored some soil freezing or thawing events because the soil at the 
measuring times was unlikely to have reached the minimum and maximum temperatures. 
AMSR-E for this study passed over the study regions at 01:30 and 13:30 when the soil 
temperatures tend to be close to the maximum and minimum values (Xie et al., 2013; Chai et al., 
2014). It means that it can capture real daily soil freezing and thawing events. Additionally, the 
results from Anandhi et al. (2013) and Wang et al. (2015) were based on data measured at 
individual sites with large uncertainties and unrepresentativeness at a regional scale because of 
their sparse and unevenly distribution compared with the satellite remote sensing data. Moreover, 
extensive human activities may also have played a role in the near-surface soil freeze-thaw cycle. 
Wang et al. (2015) reported that the observed decrease in the duration of freezing in areas with 
high rates of urban expansion was approximately higher by 42% to 52% than those of areas 
having low and medium rates of urban expansion, being consistent with our results. 
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Fig. 9 Variations in the annual air temperature during the period 2002-2011 over the SRYR. ATave, the annual 
average air temperature; ATmax, the annual maximum air temperature; ATmin, the annual minimum air temperature; 
ATspring, the average air temperature in spring; A7Tsummer, the average air temperature in summer; ATautumn, the 
average air temperature in autumn; and ATwinter, the average air temperature in winter. 


Table 1 Correlation coefficients of the duration of stages of soil freeze-thaw cycles and air temperature for the 


SRYR 


Duration of 
complete freezing 


Duration of the autumn 
freeze-thaw transition 


Duration of 
complete thawing 


Duration of the spring 


Correlation coefficient s 
freeze-thaw transition 


AT ave -0.145 0.423 0.281 0.757" 
AT max 0.065 -0.059 0.310 -0.355 
AT min -0.486 0.743" -0.450 -0.819 
AT spring 0.601* 0.236 -0.282 -0.451 
AT iina -0.259 0.562" 0.392 ~0.495 
AT autumn -0.213 0.395 0.551" 0.568" 
AT winter 0.666" 0.434 0.295 0.842" 


Notes: **, Correlation coefficient at the 0.01 level; *, Correlation coefficient at the 0.05 level. 


The variations in near-surface soil freeze-thaw cycles could significantly impact the ecosystems 
of the SRYR. The earlier thawing and the increase in duration of thawing may result in the earlier 
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onset and lengthening of vegetation growing seasons (McDonald et al., 2004; Delbart et al., 2008; 
Richardson et al., 2013). However, the increase in duration of thawing, together with climate 
warming, has a negative impact on ecosystems within water-limited regions through increasing 
water stress (Seneviratne et al., 2010; Zhang et al., 2011). Moreover, changes in soil freeze-thaw 
cycles would intensify the water and energy cycles through soil ice-water phase changes and 
increase in the exchange of heat between the atmosphere and the ground surface. Such variations 
tend to generate drier land surfaces because both water infiltration and evapotranspiration 
increase under these conditions in the SRYR (Jin et al., 2009; Li et al., 2012). A lower 
groundwater table and a decrease in the areas of lakes and wetlands have also occurred on the 
plateau, which has further accelerated the degradation of alpine grassland and desertification. 
Additionally, it has been predicted that the temperature of the Qinghai-Tibetan Plateau could 
increase by 2.2—2.6°C by 2050, with an annual increase rate of 0.044—0.052°C/a (Qin et al, 2002). 
Under the highest warming rate of 0.052°C/a, the permafrost area would decrease by 13.5% (Nan 
et al., 2005), and the mean active layer thickness would increase by 0.3—0.8 m compared with that 
in 2010 (Zhang and Wu, 2012). At the same time, the near-surface soil freeze-thaw cycle would 
severely change in the future under those climate scenarios. 


5 Conclusions 


In this study, we applied a dual-indices algorithm in the SRYR based on measured LSTs and 
vertically polarized 18.7 and 36.5 GHz AMSR-E brightness temperature data obtained through 
passive microwave remote sensing. The daytime and nighttime brightness temperature thresholds 
of 257.59 and 261.28 K used in the algorithm were confirmed for distinguishing the near-surface 
freeze-thaw status. Additionally, we divided the inferred status into three types (DFS, DTS and 
DFTS) and separated the annual freeze-thaw cycle into four stages. Most of the SRYR 
experienced freeze-thaw cycles of long duration. In the coldest month, the DFS occurred 
primarily in area I, DFTS occurred in areas III and IV, and DTS appeared in area II. The areas 
occupied by DFS and DFTS during the period 2002-2011 decreased by 0.08% and 0.25%, 
respectively, and the DTS area increased by 0.36%. According to the east-west profile line at 
34.5°N and north-south profile line at 99.8°E, the onset of thawing advanced, the onset of 
freezing became delayed, and incomplete freezing occurred on the Zoige plateau in the 
southeastern part of the SRYR. During the period 2002-2011, the dates of the onset of thawing 
and complete thawing advanced by 3.1 (+1.4) and 2.46 (+1.4) days, respectively; the dates of the 
onset of freezing and complete freezing delayed by 0.9 (+1.4) and 1.6 (+1.1) days, respectively. 
The duration of complete thawing increased by 0.72 (+0.21) day/a and the duration of complete 
freezing decreased by 0.52 (+0.26) day/a. Most severe changes occurred in areas II and IV. The 
durations of the spring and autumn freeze-thaw transitions increased over most of the central 
region, whereas the durations decreased in other regions. 

The results of this study demonstrated that the increases in ATmin and AT winter were the primary 
drivers for earlier thawing and later freezing and also for an increase in the duration of thawing 
and a decrease in the duration of freezing. Given the fragile and sensitive ecosystems that occur 
across the SRYR, these variations in the near-surface freeze-thaw cycle affect the water and 
energy cycles in the soil and cause further acceleration of the degradation of grasslands. 
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